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';..  •  i    I  INTRODUCTION 

The  Weibvill  distribution  was  originally  derived  in  a  paper  by  Fisher 
and  Tippett  (2).  However,  as  the  name  implies,  it  is  referred  to  as  the 
Weibull  distribution  because  of  a  paper  by  Weibull  (8J  some  10  yeejrs  later. 
It  is  a  three  peirameter  distribution  but  without  undue  loss  of  generality, 
the  location  pareimeter  c  can  be  set  equeJ.  to  zero.  The  two-parameter  caise 
will  be  considered  in  this  report. 

The  Weibull  distribution  plays  an  important  role  in  describing  many 
natvtreil  phenomenon.  Some  examples  given  in  Weibxill  (9]  are:   (l)  Yield 
strength  of  steel,  (2)  Fiber  strength  of  Indian  cotton,  (3)  Statvires  of 
adult  males,  born  in  the  British  Isles.  Dubay  (l)  describes  how  the 
Weibull  can  be  used  to  describe  the  deterioration  pattern  of  certain 
inexpensive  industrial  products.  In  his  examples,  time  to  deterioration  is 
the  random  variable  which  obeys  a  Weibull  probability  law. 

One  of  the  difficiilties  encountered  in  working  with  the  Weibiill  dis- 
tribution involves  estimating  its  parameters.  The  purpose  of  this  report 

is  to  show  four  methods  of  estimating  the  shape  and  scale  pareuneters  of 
the  Weibull  distribution;  namely,  the  method  of  moments,  majcimum  likelihood 
method,  variance  method,  and  Rand  Corporation's  method.  These  will  be 
derived  in  some  detail  and  properties  of  these  estimators  will  be  included, 
when  possible.  '  :  -   • 


I      ■.   II  CHARACTERISTICS  OF  THE  WEIBULL  DISTRIBUTION 

A  random  variable  X  is  said  to  distributed  as  the  Weibull  distri- 
bution if  its  density  is: 

f(x;a,b)  =  abx^"-*-  exp{-b(x)^}  ;  X  >  0      ^^       (l) 

•     ■     •  =0  X  <  0   . 

This  is  a  two  parameter  feunily  of  distributions,  the  parameters  being  a 
and  b.  Both  the  shape  parameter  a  and  the  scale  parameter  b  must  be 
greater  than  zero.  Figvires  1,  2,  and  3  show  graphs  of  the  Weibull  distri- 
bution for  various  values  of  a  emd  b.  Figure  1  shows  the  distribution 
to  be  exponential,  when  both  a  and  b  are  set  equal  to  1.  Figure  3  shows 
the  distribution  to  be  "almost"  bell-shaped  when  a  =  3.26  and  b  =  O.TO. 

The  mean,  variance  and  median  will  be  derived  for  the  density  given  in 
(1).  To  derive  these  qviantities,  one  needs  to  find  /  xf(x)  dx,  /  x  f(x;  dx 
euid  the  value  m  such  that  jr^f(x)  dx  =  1/2  respectively.  By  vising  the 
trems formation  y  =  bx  ,  ' 


^"xf(x)  dx  =  /"(y/b)-'-^*  •  exp(-y)  dy 


-1/a  .  »  (y)^/^ 


=  rd+l/a)  .  b-  '   .  /  fLL-^   .  exp(-y)  dy 

=  r(i+i/a)  •  b"-'-/* 


•  1. 

because  the  above  integrand  is  just  a  gamma  density  function.     Hence  the 
meeua  of  the  Weibvill  distribution  is  given  by: 


i . 


y  =  h-^^^  .   r(l+l/a)   .  (2) 
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By  using  the  transformation  y  =  bx  again 


^Vf(x)  dx  =  /"(y/b)^''*  •  exp(-y)  dy 

^     ,  =  r(l+2/a)  •  b"^^^  •  f  y(/^2/&)    '   ®^^-y^  ^ 

=  b"^/*  .  r(l+2/a)  ;  1  (3) 

because  the  integrand  is  just  a  gamma  density  function  again.  The  variance 
of  a  random  variable  is  E(X^)  -  E^(X).  Thus,  the  variance  of  the  Weibull 
distribution  is  given  by: 

;  o2  =  b"^/*  {r(l+2/a)  -  r2(l+l/a)}  .  (U) 

Finally,  it  can  be  shown  that   . 

:  ^-l/a(ln  2)^/* 

r  f(x)  dx  =  1/2 

(by  xising  the  transformation  y  =  bx*  again).     Hence  the  median,  m,   is  given 

by:  ■    ■     .     ,  ■"'■  ; 

t..:    .     ■;  '  m  =  b-^/^ln  2)^/^     . 

In  a  seeirch  for  "good"  estimators  for  a  and  b  in  (l),  one  should 
consider  the  property  of  sufficiency.  Sufficiency  is  an  especially  important 
property  because  if  sufficient  estimators  exist  then  one  can  proceed  to  find 
minimum  varianced  unbiased  estimators.  In  general.,  an  unbiased  estimator 
baised  on  a  sufficient  statistic  will  have  smetller  variance  than  one  which 
is  not  based  on  a  sufficient  statistic  (see  Mood-Graybill  [Tj »  P  176). 


Unfortunately,  a  set  of  Joint-sufficient  statistics  does  not  exist. 
To  show  this,  consider  that  in  a  random  sample  of  size  n  (n  >^  2),  the 
likelihood  is  given  by:         . 

;  '   ;  ':  ■     L  =  bV  n  x^"^  exp(-b  ^  x^)  (5) 

i=l  i=l  ■'■ 

In  order  to  obtain  a  set  of  joint-sufficient  statistics  for  a  and  b, 
say  6^^  euad  e_,  the  likelihood  must  be  expressible  as: 

L  =  k(x^,...,x^)  g(a;b,9^,e2) 
vhere  k  does  not  contain  a  or  b.  Consider  the  term  in  (5), 


exp(-bx  )  exp(-bx  )  .  .  .  exp(-bx  )  . 


Clearly,  it  is  impossible  to  separate  the  x. 's  from  a.  Thus  joint- 
sufficient  statistics  for  a  and  b  do  not  exist.  ; 
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•;..-  Ill  METHOD  OF  MOMENTS 

This  method  of  estimation  simply  equates  the  first  two  saurple  moments 
vith  the  first  two  population  moments.  The  first  two  population  moments  as 
derived  in  II  £ire: 


^ 


=  E(x)  =  ■b"^''^  r(l+l/a). 


yg  =  E(x^)  =  b-2/*  r(l+2/a)  . 


The  seuaple  moments  are  given  hy: 


n 
m,   =     J]     x./n  =  X,  (6) 

■^       1=1     ^ 


,  .  .  .  m^  ""  I  ^i/"  •  ^'^^ 

Now  proceeding  to  derive  the  moment  estimators  by  the  usual,  method,  let 


f 


Pg  =  m^. 


Then 


x=b-^/*  r(l+l/a),  (8) 


n 

I 
i=l 


I   x^/n  =  b-2/^P(l+2/a)  .  (9) 


The  quantity  b  can  be  eliminated  by  dividing  equation  (9)  by  equation  (8) 
squared  which  yields: 


I  X  /n 
1=1      _  r(l->-2/a) 


-2 

X 


r'^d+i/a) 


(10) 


Clearly,  (lO)  ceuinot  be  solved  explicitly  for  a,  but  with  the  aid  of 
tables  of  the  Gamma  function  one  can  determine  the  value  of  a  which  sat- 
isfies (10).  This  is  the  method  of  moments  estimator  for  a. 

After  a  has  been  determined,  equation  (8)  csui  be  solved  for  b. 
The  solution  is  the  method  of  moments  estimator  for  b  and  is  given  by: 


b  =  {(1/x)  •  r(i+i/a)r  . 


(11) 
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IV     METHOD  OF  MAXIMUM  LIKELIHOOD 


The  likelihood  fvinction  for  a  reuidom  sample  x,  ,x^,....x     from  a  Weibvill 

1  2     n 

distribution  can  be  written  as 

n     ,  n 

a-1 


L  =  e!\^     n  x^~  •  exp{-b(  ^   xf)} 
i=l  ^  i=l  ^ 


or  upon  taking  the  naperian  logarithm  of  L,  one  obtains, 


«  n  n 

•  L  =  In  L  =  n  ln(b)  +  n  ln(a)  +  (a-l)  I     ln(x. )  -  b  I     xt  . 

i=l  i=l 

« 

To  get  the  maximiom  likelihood  estimators,  L  is  differentiated  with 

respect  to  a  and  then  with  respect  to  b.  The  resulting  equations  are 
set  equal  to  zero  and  solved  for  a  and  b.  Thus 


~-=^^     I     ln(x.)  -b  I     x^ln(x.)  . 
da   a   . ''t     1      .'',1     1  ' 
1=1  1=1 


dL    n 


db  "  b  "  ^^,  ^i  ' 


which  yield  the  following  equations  in  a  and  b. 


n 


^+  I     ln(x.)  -b  I     x^ln(x.)  =0, 
a   i=l    ■"     i=l  ^     ^  . 


(12) 


--  I  x^=0  .   ■  (13) 

b   1=1  ^ 


11 


However,  these  equations  cannot  "be  solved  explicitly  for  a  and  b, 
(see  Menon  (6)).  In  fact  (12)  cannot  be  solved  for  a  even  if  b  were 

X     X 

known.  Intuitively,  the  reason  is  that  equations  of  the  form  d  +  e  =  k 
cannot  be  solved  explicitly  for  x  where  d,  e,  and  k  are  constants. 
Thus  an  iterative  scheme  is  sought  which  will  yield  the  values  of  a  and 
b  in  equations  (12)  and  (13). 

The  iterative  technique  developed  by  Hardy  (3)  is  to  write  equations 
(12)  and  (13)  as:        • 


n 
a 


J+1      ^       ^^  ^ 

.   b,  I     x  J  ln(x  )  -  I     ln(x  ) 
J  i=l  ^     ^    i=l     ^ 


> 


(11.) 


',  °   °  „   •  (15) 


Now  guess  an  initieQ.  value  of  a,  call  it  a  .  This  is  put  in  equation  (15 ) 

to  yield  a  value  of  b  .  These  two  values  of  a   aj:id  b   are  substituted 

o  00 

in  (li*)  to  yield  a  value  of  a,  .     This  process  is  repeated  until     a.     and    b. 

converge.     These  convergent  values,     a.     emd    b       are  the  meiximum  likelihood 

J       J 

estimators  for  a  and  b  respectively.  Lloyd  and  Lipow  [k]  claim  that 
this  method  will  converge  if  a  reasonable  estimate  of  a  is  made.  The 
moment  estimator  for  a  given  in  (lO)  can  be  used  ais  the  initial  value  a  . 


12 


V  RAND  CORPORATION'S  METHOD 

This  method  of  estimation  is  due  to  Wilson  [lo) .  A  bayesian  approach 
is  used;  namely,  it  is  assumed  that  the  scale  parameter  b  has  a  prior 
gamma  distribution.  Hence  one  can  obtain  the  joint  density  of  x^,...,x^ 
and  b  and  then  integrate  this  Joint  density  with  respect  to  b  to  obtain 
a  marginal  likelihood  function  which  contains  only  the  shape  parameter  a. 
Finally  asymptotic  properties  of  this  likelihood  function  yield  an  estimator 

for  a. 

In  a  random  sample  of  size  n,  the  joint  density  of  the  x's  given  b 

is:  ;■.,  '■■^'y- 

n   n(a-l) 


h(x^,...,xjb)  =.a\V^^-^^  exp(-bnr(a)^) 


where 


g  =  n  X. 
■  i=l  ■ 


1/n 


and 


n 


i=l 


r(a)  = 
Now  ass\ime  the  density  for  b  to  be  given  by: 


1/a 


g(b)  =« 


,.^N' 


i&yb^'-^  exp(-bT')  ;  b  >  0 


;  b  <  0 


Since 
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f(b,Xj^^,    ...»  xj  »  h(x^,    ....  x^|b)   g(b). 


/"f(b,x,,    ...»  X  )db  a  /"f  db     will  contain  only     a,  T' , 


and  N'    . 

Now 


N' 
/••*jv        f"  ^^  n(a-l)  /  ,       /    \a\    (T' )         ^N'-l  ,  vm^^v 

^  fdb  =  ^  al)  g  '  exp(-bnr(a)    )    ^p/^,k     b  exp(-bT')db, 


n  n(a-l)    (T' )^'    ,«  n+N'-l         ,  ,  ,^,^     ,    ^a^, 
ig  r(H')   o  exp{-b(T'+nr(a)    )} 


Therefore 


-»   -,  n  n(a-l)/„,xN'    r(N'+n)       1 ,,^v 

^fdb  =  ag  (T-)       -TTFT  •    ,^.,     I    ,a.N'^n  ^1^) 


t    I 


[T'+nr(a)**] 
since  the  integrand  can  be  put  in  the  form 


;     '  N'+n 

fT'+nr(a)^)  ,n+N'-l  ,  v/mt^     t    sQ-w 

r(H'-»ni    *   exp{-b(T'+nr(a)    )} 

which  is  a  gamma  density  function.  Hence  (l6)  is  the  marginal  likelihood 
function  for  a,  call  it  L(a).  From  (l6),  L(a)  is  proportional  to 

a°(g)'^^  [l/(T'+nr(a)^))  "''.-.  (17) 

From  (17),  one  can  obtain  the  maximiom  likelihood  estimator  directly; 
however,  it  is  a  burdensome  computation.  Therefore  asymptotic  properties 
of  L(a)  are  exploited  to  eurrive  at  an  estimator  for  a. 


11* 


A  Theorem  from  Mann  and  Wald  (5)  states  that  if  a  sequence  of  random 
numbers  converge  in  the  probability  sense  to  c,  then  g(x^)  converges  to 
g(c)  for  any  continuous  function  g.  It  is  used  to  prove  the  following: 
Lemma  1:  . 

p  lim  r(a)  =  ,  ^^]\ 

n-«  '  '   r(l+l/a)  . 

where  p  lim  defines  a  limit  in  the  probability  sense. 

Proof:  V 

From  (2).  E(x)  =  b"^/^  r(l+l/a).  By  applying  the  Weak  Law  of  Large 
Numbers  to  E(X)  the  following  is  obtained: 

'   .   P  lis  X  =  b-1/^  r(l^l/a)  ,  '         (18) 

or  since  X  =  r(l), 


p  lim  r(l)  =  b-1/*  r(l+l/a)  . 


^^   xva.,  -  u     iu+J./a;  .  (15) 


Clearly,  the  distribution  of  W  =  X*  is  exponential  with  mean  given  by 

1/b.  Again  by  applying  the  Weak  Law  of  Large  Numbers  to  E(W)  the  following 


is  obtadned: 


•:   :;  ;  ■    ■       p  limW=  1/b  . 

Now  consider  the  function  W^^^  (which  is  just  r(a)) 
Hence,  by  the  Mann  Wald  Hieorem  the  above  equation  becomes: 

,  1;   •        Pli2r(a)  =  (l/b)^/^  .  (20) 

Combining  (19)  and  (20),  it  follows  that 
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PjiS^^^^  "  Tilil/a) 


This  completes  the  proof  of  Lemma  1. 

To  obtain  an  asymptotic  expression  of  L(a)  in  (17),  one  can  use 
Lemma  1  and  the  Mann  Wald  Theorem  to  write  (l?)  as: 


L(a) 


n  na 
a  g 


T'+n 


r(l) 

r(i+i/a) 


N'+n 


+  Error  term.   (21) 


Now  for  practical  reasons  the  major  result  is  the  following. 


Lemma  2: 


p  lim  f-  In  (L(a))  +  In  n]  =  ln(a)  +  a  In 


n-**»  ^n 


g  •  rd-n/a) 


<1) 


(22) 


Proof: 

The  proof  follows  immediately  if  one  takes  the  logarithm  of  both 
Bides  of  (22).  This  operation  yields  the  following: 


-  In  (L(a))  +  ln(n)  =  ln(a)  +  a  In  g 
n 


(1  +  — )  in 
n 


N' 


T'/n  + 


■iiy 


r(i+i/a)' 


ln(n)  +  Error  term  . 

n 


Now 


take  p  lim  of  both  sides  to  get 
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p  lim  [-  ln(L(a))  +  ln(n)]   =  ln(a) 


+  a  In 


K  .  rd-n/a) 


The  maximum  asjinptotic-marginal-likelihood  estimate  of  the  shape 
parameter  is  defined  £is  the  value  of  a  maximizing  the  likelihood  function 
(22).  Note  that  this  estimator  must  be  distinguished  from  the  asymptotic 
form  of  the  majcimum  margined,  likelihood  estimator  that  can  he  derived  explicitly 
from  (16). 

If  the  derivative  of  the  right  hand  member  of  (22)  with  respect  to 
a  is  set  equal  to  zero  the  following  implicit  relation  for  a  results: 

(i|;(l/a)  -  1)   a  -  In  (r(l+l/a))  =  ln(g/r(l))  , 

where  ';    . 

t|»(x)  =  r'(x+i)/r(x+i)  . 

Values  of  g/r(l)  are  tabulated  in  Table  1  for  selected  values  of 

n  "'"■^' 

a  [recstll  r(l)  =  I     x. /n]  .  This  table  is  a  reprint  from  Wilson  fio]  . 
.  i=l  ^ ■ 

Hence  to  obtain  a  one  must  calculate  g/x  and  then  go  to  the  table  to        ;, 
find  the  value  of  a. 


Table  1 


MAXIMUM  ASYMPTOTJC-MARGINAL-LIKELIHOOD  ESTIMATES  FOR  THE  SHAPE 
PARAMETER  OF  A  WEIBULL  DISTRIBUTION 
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g/x 


g/x 


.11*59 

.0500 

.1+285 

.1496 

.0526 

.5615 

.1536 

.        .0556 

.5825 

.1580  •. 

.0588 

.6060 

..;,.         .1627    •'               > 

.0625 

.6321+ 

.1680     . 

.0667 

.6625 

1      .1737  •.     . 

.071I4 

.6970 

.1801  . 

.      .0769 

.7372 

■:.        .1873      ■ 

.0833 

.  .781+6 

4       .1953.  . 

*0909 

.8117 

.  201+6 

:•          .1000 

.81+16 

.2152 

.1111 

.8961+ 

.2278  ■ 

.1250 

.9116 

.21+28;    •• 

.11+29 

.9276 

.2612  :-\  - 

.1667 

.9531 

.281+5    ■ 

.2000    • 

.9620 

.3155 

.2500 

.9805 

.3591* 

.3333 

.9901 

.5000 

1.0000 

1.1111 
1.2500 
1.1+286 
1.6667 

2.0000 
2.5000 
3.3333 
l+.OOOO 
5.0000 
8.3333 

10.0000 
12.5000 
20.0000 
25.0000 
50.0000 
100.0000 
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VI  VARIMCE  METHOD 


In  considering  this  method  due  to  Menon  (6),  the  Weibull  density  is 
written  in  a  slightly  different  form,  namely: 

(a/b')(x/b •)*"■'•  exp{-(x/b')*};  x  >  0 


f(x)  =  < 


(23) 
;  X  <  0  . 


In  this  form  of  the  Weibiill  distribution  a  and  b*  are  the  shape  and 
scale  parameter  respectively.  The  relation  of  (l)  to  {2h)   is  that 
b  =  d/b')^.   In  this  case  a  different  shape  parameter  is  being  estimated. 
The  method  employed  is  to  find  an  estimate  of  (l/a),  say  d  =  (l/a) 

and  then  set  a  =  (l/d).  Similarly,  an  estimate  of  In  b'.  In  b'  is  obtained 
and  then  b'  is  set  equal  to  exp(ln  b*).  Of  course,  it  is  known  that  a 
and  b'  do  not  necessarily  possess  the  same  statistical  properties  that 
d  and  In  b'  possess.  Nevertheless,  the  justification  of  b'  and  a  as 

A  A 

estimators  will  be  in  the  "good  qualities"  of  In  b'  and  d. 

,   In  order  to  derive  the  estimator  a,  consider  the  random  variable 
W  =  (X/b')^.  Since  dw  =  (a/b' )(x/b' )^""^  dx,  clearly  the  distribution  of 


w  is: 


;(v)  = 


e"^    w  >  0 


w  <  0 


Hence  W  is  distributed  independently  of  a  and  b'.  In  particular 
In  W  is  independent  of  a  and  b'.  By  definition,  the  variance  of  In  W 

is  given  by  E(ln  W)  -  E  (in  W) .  Menon  (6)  has  calculated  the  value  of 
V6Lr(ln  W)  to  be  11^/6.  Also  since  W  =  (X/b')  ,  the  following  is  obtained: 
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Var  (In  W)  »  Var  {ln(X/b»)*} 


Var  {a(ln  X  -  In  b')) 


=  a  Var(ln  X)  . 
Thus  by  combining  the  two  expressions  for  Var  (in  W),  the  following  results! 


1/a  =  {Var  (in  X)  •  6/n2}^/^ 


(24) 


Now  set  d  =  1/a  emd  the  estimator  d,  of  d  (also  of  l/a)  is  obtained  by- 
putting  the  sample  variance  of  In  X  in  (2i+),  or: 


d  = 


/. 


6/n2 


.2   i=l 


n 

I     InX.)' 


I   (In  X.)'^  -i::i— 
i=l 


n-1 


1/2 


As  stated  previously,  the  estimator  of  a  denoted  a  is: 

a  =  1/d  . 

Menon  [6j ,  reports  the  following  statisticaJ.  properties  of  d.  The 

quantity  d  is  not  an  unbiased  estimator,  however,  its  expected  value  and 
variance  is  given  by: 

■ 'E(d)  =  d  +  0(l/n) 


2 

Var  (d)  =  1.1  —  +  d^  0(l/n^)  . 


For  large  n,  d  is  asymptotically  normal  with  meeui  and  variance  given 
above.  Its  asymptotic  efficiency  is  55^. 

One  is  really  interested  in  the  properties  of  a.  Not  much  is  known 
about  the  statistical  properties  of  a.  By  Menon  (6]  ,  a  is 
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•.       '■     ,  ■  2 

asymptotically  normal  with  mean  and  variance  given  by  a  and  1.1  — 

respectively. 

The  estimator  b'  can  be  obtained  after  the  value  of  d  from  above 
is  known  by  using  the  following  argument.  Consider: 

.',:    E(ln(X))  -  /"in  x  (^)(^)    exp{-(^)  }  dx 

'a  -,   a-1 

Let  y  =  (^)   ,  then  dy  =  (fT)(^)     and  E(ln(X))  becomes: 


E(ln(X))  =  rin(b'  •  y^/*)  e'^  dy 
o 


=  /"(inb-)  e-y  +  rin(y^/^)  •  e'^dy 


=  in  b '  +  -  /"in  y  •  e"^  dy  . 
a  o    "' 

Recall  that  d  =  1/a  and  since  /*(ln  y)  e~'^  dy  =  -.5772  the  above  equation 
becomes:  •' 


E(ln(X))  =  In  b'  +  d(-.5772)  . 

n 
Now,  in  the  above  equation  one  can  estimate  E(ln(X))  by  1   Inx./n,  In  b' 

i=l   ^ 

by  In  b'  and  d  by  d  to  get  the  following  estimator  for  In  b': 


«     n 
In  b'  =  j;  In  X.  +  d  (.5772), 
i=l    ^ 

or  upon  substituting  the  value  of  d  obtained  previoiisly  the  following 

estimator  for  In  b  is  obtained: 


*     -  .r.'r 


n 


iix  b'  «  I  la  X Jn  +   .5772 

•  ■.  i=l     . 


n 


(6/n2)  {(  I   (lnx.)2  --^^ 
i=l    ^ 


n      i 
illnx.) 


n 


)/n-l} 
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1/2 


Now  it  seems  reasonable  to  define  the  estimator  for  b'  to  be 

.  b'  =  expdn  b')  . 
Menon  [6]  reports  the  following  statistical  properties  of  In  b' 

E(ln  b' )  =  In  b  +  0(l/n) 


Var  (In  b')  =  1.2  d^/n  +  d^  0(n"^/^)  . 

The  quantity  In  b'  is  asymptotical  normal  with  mean  and  variance  given 
above.  Again,  little  is  known  about  the  statistical  properties  of  the 
estimator  for  b',  exp(ln  b').  But  by  Menon  (6) —  b'  is  asymptotical  normal 

with  mean  and  variance  given  by  b'  and  1.2  dHa  /n  respectively. 
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VII  CONCLUSION 

Four  methods  of  estimating  the  shape  and  scale  parameters  of  the 
Weibull  distribution  were  presented.  It  was  shown  at  the  onset  that  suf- 
ficient estimators  do  not  exist.  The  estimators  derived  by  the  variance 
method  were  shown  to  be  asymptotically  normal.  With  this  exception,  prop- 
erties of  the  proposed  estimators  are  a  subject  of  further  investigation. 

For  further  study  it  is  recommended  that  the  properties  of  these  esti- 
mators such  as  bias,  variance,  consistency,  and  efficiency  be  studied, 
perhaps  by  Monte  Ceurlo  methods. 
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A  random  variable  X  is  said  to  be  distributed  as  the  Weibull  distri- 
bution if  its  density  is  given  by: 

•   f(x;a,b)  =  abx*"^exp[-b(x)^)   ;  x  >,  0 

=  0  ;  X  <  0  . 

This  is  a  two  parameter  family  of  distributions  in  which  the  shape  parameter 
a  and  scsLLe  parameter  b  must  both  be  greater  theui  zero.   In  its  most 
genereil  form  the  Weibull  distribution  contains  a  location  parameter,  but 
without  loss  of  generality  the  location  parameter  is  assumed  to  be  zero. 

Four  methods  of  estimating  the  parameters  of  the  Weibull  distribution 
6Lre  presented:  the  method  of  moments,  maximum  likelihood  method.  Rand 
Corporation  method,  and  the  variance  method. 

The  method  of  moments  estimators  are  derived  by  equating  the  first  two 
popvilation  moments  about  the  origin  with  the  first  two  sample  moments  about 
the  origin.  This  procedure  yields  an  implicit  equation  for  a  which  may 
be  solved  when  the  sample  values  are  known.  When  a  is  determined  the    / 
explicit  form  of  b  is  given. 

The  maximum  likelihood  estimators  are  derived  by  differentiating  the 
likelihood  function  with  respect  to  a  and  b.  However  when  the  resulting 
equations  are  set  equal  to  zero,  explicit  expressions  for  a  and  b  do 
not  exist.  But  by  iising  iterative  techniques,  solutions  for  a  and  b 
can  be  found. 

Rand  Corporation  develops  an  estimator  for  the  shape  parajneter  a 
by  vising  a  bayesietn  technique.   It  is  assumed  that  the  scaile  parameter 

b  has  a  prior  gamma  density.  Then  the  Joint  density  of  x.,....,x  and  b 

1      n 

is  obtained  and  integrated  with  respect  to  b  to  obtain  a  marginal 


likelihood  function  which  contains  only  the  shape  parajneter  a.  Asymptotic 
properties  of  this  likelihood  fvinction  yield  an  implicit  equation  in  a. 

A 

However  with  the  aid  of  Table  1,  a  can  be  determined  after  the  sample  data 
is  known.    • 

A  A 

The  variance  method  first  obtains  an  estimate  d  of  l/a  and  then  the 
estimate  of  a  is  given  by  a  =  1/d  .  Similarly  an  estimate  In  b  is  derived 

A  A  A 

and  then  the  estimator  of  b  is  given  by  b  =  exp  (in  b).  Both  a  and 

A 

d  are  asymptotically  normal. 


